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A scattering method is used to calculate the surface band structure of AI(lll) from 8.6 cV below 
the Fermi level to 9 cV above it. This method has rarely been implemented previously. The complete 
complex bulk and surface band structure is also shown for key values of the crystal surface-parallel 
momentum. The surface region is treated semi-empirically from a small number of experimentally 
determined quantities. We find agreement over the whole surface Brillouin zone for surface states 
compared with recent theoretical ab initio results and experiment. Differences with respect to 
surface resonances detected by recent theoretical methods are discussed. In particular a class of 
surface resonances is detected here that may be surface-layer localised and results from very wide 
partial surface-projected energy gaps where real lines in the complex bulk band structure extend 
from band minima to infinity. Other metfiods do not appear to have detected these resonances. They 
may have an effect on photoemission intensity profiles and may account for previous unexplained 
experiment features considered to be "a kind of surface enhancement at partial gap edges" . This 
scattering method provides a single picture for the formation and classification of surface states and 
resonances. 
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I. INTRODUCTION 

We use the scattering matrix method of McRae-'^ to cal- 
culate the surface band structure E{ki\) for semi-infinite 
Al(lll) to 9 eV above the Fermi energy for M - K - 
r - M of the surface Brillouin zone and show the com- 
plete near-surface energy band structure E{k\\,ki_) for 
symmetry points F, M, K and some points along the di- 
rection r - Af(E). These complete E{k\\ , k±) band struc- 
tures show bulk states with propagating and attenuating 
wave functions as well as all surface states and resonances 
(SS/SRs) that are needed in the interpretation of photoe- 
mission and other surface sensitive phenomena. Reliable 
methods for clean metal siirfaces are needed before calcu- 
lations for more complex systems of foreign-atom layers 
on metal substrates can be attempted. 

There have been a number of calculations of band 
structure for this surface by various methods^"^° and re- 
cent methods*"^*^ including the present one agree on the 
occurrence of surface states. All the SS/SRs that have 
been found experimentally to-date are found in the re- 
cent calculations and present one. However even the re- 
cent ah initio self-consistent calculations with free-form 
surface potentials in Refs. 8-10 do not give consistent re- 
sults for SRs and this does not inspire confidence in the 
reliability of any of the calculations for this case. The 
surface-embedded Green-function (SEGF)^ method iden- 
tifies SS/SRs from the calculation of the layer-density- 
of-states (LDOS). The repeated-slab density-functional 
(DF) method^ identifies the SS/SRs from maxima near 
the surface and decay into the bulk of the planar charge 
density as a function of the coordinate perpendicular to 
the surface and also from the LDOS. The assembly of 
boundary-controlled monolayers (ABCM)'^^ method re- 
lies on matching attenuating crystal wave functions to 



vacuum wave functions. The present scattering method 
identifies SS/SRs from a condition for confinement of the 
electron in the surface region due to sustaiiicnl multiple 
scattering in that region. Our result agrees very closely 
with that of ABCM for a particular class of surface res- 
onances. In the present scattering method we also find 
another class of surface resonances that would not be de- 
tected in the ABCM with the search method used in that 
particular case. Some of these "extra" resonances but 
not all, appear in the slab DF calculation. If reasons can 
be found for the different results between methods then 
the most reliable methods can be identified. Agreement 
between calculated surface band structures from differ- 
ent theoretical approaches would provide confidence in 
the results and aid in the analysis of measured features 
found in inverse and direct photoemission. 

From our calculation we suggest that the resonances 
that are missing in the work of Ref. 10 occur in surface- 
projected bulk-band gaps or partial gaps that are very 
wide in energy and contain "real lines" with real energy 
in the complex band structure where the minima of bulk 
bands have no lower-energy band maxima to which to 
join and extend to infinity. This type of resonance is of- 
ten very surface-layer localized and may have been missed 
in the search method of Ref. 10. Additionally it is sug- 
gested that unexplained features in experimental angle- 
resolved photoemission spectroscopy (ARPES) that have 
previously been termed "partial bulk gap-edge features" 
arising from "a kind of surface enhancement at partial 
gap edges" are surface resonances of this type as these 
resonances often occur at energies just outside other nar- 
rower absolute or partial projected bulk gap-edges that 
lie within the energy range of the very wide gap. Such 
resonances are just as easy to detect using the present 
method as are SSs. It is suggested that this property 
gives the method some advantage over other methods. 
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The above recent calculations are ab initio and deter- 
mine the surface potential self-consistently from a free- 
form shape. In most cases they use a local-density ap- 
proximation for the exchange-correlation potential and 
do not produce the image potential in the vacuum 
that varies as the inverse of distance from the surface. 
The present method determines the non-spherically sym- 
metric surface potential semi-empirically and other key 
properties from experiment. The geometric structure 
of the surface layers is determined from the result of 
previous low-energy electron diffraction (LEED) anal- 
ysis. The spatial variation of the valence charge den- 
sity at the crystal-vacuum interface and hence the sur- 
face barrier potential is found from previously deter- 
mined experimental energies of image barrier and other 
SS/SRs at r and the experimental work function. The 
spherically-symmetric bulk ion-core potentials and con- 
stant interstitial potential are determined from previous 
self-consistent bulk-band structure calculations that have 
shown agreement with experimental bulk states. Experi- 
mental surface core-level binding energy shifts for Al(lll) 
are negligible^^ and final-state relaxation contributions 
are small for simple metals. This leads to choosing 
the surface-layer constant potential and the surface ion- 
core potential the same as the bulk in this case. The 
non-spherically symmetric surface potential is then ef- 
fectively an overlap of the surface barrier potential and 
spherically- symmetric surface ion-core potential. In this 
scattering formulation ion-core potentials are represented 
as point sources once their scattering phase shifts are cal- 
culated and their radius is no longer of relevance. It is to 
be noted that there are no discontinuities in the constant 
potential or its derivative in the surface region in these 
calculations that might give rise to spurious non-physical 
features. 

While this surface potential may seem too crude an 

approximation in comparison with the ab initio self- 
consistent surface potentials of the above recent calcula- 
tions it adequately produces all SS/SRs found from those 
methods in this case. It can also be used at energies above 
the vacuum level where inelastic effects become signifi- 
cant. Its applicability to erases wlic^re surface core levels 
shifts are expected and for overlayers of foreign atoms at 
the surface is to be assessed in the future. It may have 
use for systems that are too complicated for the ab initio 
methods. 



II. METHOD OF CALCULATION 

A scattering method is used to determine the electronic 
properties. The semi-crystal is divided into atomic layers 
parallel to the surface. Between these layers are infinitesi- 
mally thin 2-dimensional (2D) periodic layers of constant 
potential where the wave function may be expanded in 
terms of a sufficient number of forward and backward 
travelling plane waves. A plane wave inside the crystal 
scatters from the top atomic layer with energy Ei rel- 



ative to the constant potential energy between ion-core 
potential wells in the layer, Ec^i and with surface-parallel 
momentum or wave vector fey. In the general case the 
bulk layers may have a different constant potential en- 
ergy, Ec and the energy with respect to Ec is E. In 
the present case for Al(lll) we have Ec^i = Ec that will 
be justified later. The wave function of a layer is ex- 
panded into spherical waves and matched to plane waves 
(propagating and evanescent or attenuating in space) in 
the regions of constant potential between the layers us- 
ing the layer-Korringa-Kohn-Rostocker (LKKR) Green 
function method of Kambe.-'^^ Partial- wave phase shifts 
with angular momentum I describe the ion-core scatter- 
ing. It should be noted that rapid convergence is ob- 
tained in this calculation for these energies because sums 
over both direct and reciprocal net vectors are performed 
in the method of Kambe. The scattering properties of a 
layer are described by a layer transfer matrix that con- 
tains complex amplitude reflection and transmission co- 
efficients. 

The method of McRae^^ is used to combine layer scat- 
tering properties to obtain the scattering properties of 
the semi-infinite crystal and hence the electronic prop- 
erties. For bulk layers that become the substrate of the 
semi-crystal, the layer transfer matrix is combined with 
the bulk geometrical structure to form the reduced trans- 
fer matrix R of the crystal substrate. From the eigen- 
vectors of R the matrix of complex amplitude reflection 
coefficients for the whole semi-infimte crystal substrate 
Ml is obtained. Further transfer matrices for special 
surface atomic layers may be combined with Mi to form 
the matrix M that contains the complex amplitude re- 
flection coefficients for scattering from all atomic layers of 
the semi-infinite crystal. The scattering properties of the 
top-most surface barrier layer of potential are contained 
in the matrix S. The sub-matrix S^^ of S contains the 
complex amplitude reflection coefficients that describe in- 
ternal reflection at the surface barrier from inside the 
crystal surface. This is found here for a one-dimensional 
barrier potential perpendicular to the surface of arbitrary 
shape by numerical solution of the wave equation. 

In the thin layer of constant potential between the 
top atomic layer and the surface barrier are forward and 
backward travelling plane waves as described above that 
scatter between the barrier and crystal. In the simpli- 
fied casc"'^ of a single plane wave scattering in a crystal 
surface-projected band gap, the modulus of the crystal 
and barrier reflection coefficients is 1 indicating total re- 
flection. Multiple scattering between crystal and barrier 
can occur and a standing wave can form for certain val- 
ues of (£^, fe||). This corresponds to the electron being 
confined in this region in a surface state. The condition 
for the formation of a standing wave is 

(t>b + (t>c = 27rn with n integer (1) 

where (jib and (t)c are the phase change on scattering at 

the barrier and crystal respectively. In terms of complex 
amplitude reflection coefficients of barrier pb, and crystal 
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Pc, this condition is 

1 - pbPc = 0. (2) 

where pi, = exp(z<^6) and pc = exp{icp(.). For {E, fey) near 
surface-projected bulk-gap edges <pc varies rapidly and in 
the barrier image tail (j)b changes rapidly allowing Eqs. (1) 
and (2) to be satisfied often in these cases. Surface res- 
onances may occur when the modulus of the crystal re- 
flection coefficient is less than 1 giving weaker reflection, 
and scattering does not occur in a surface-projected band 
gap. In this case Eq. (2) becomes 

|1 — PftPcl^ minimum. (3) 

In the case of real crystals more than one plane wave 
must be considered and a matrix representation is re- 
quired where the reflection coefficients of barrier and 
crystal are contained in the matrices S^^ and M respec- 
tively. McRae has determined^ that in this case one or 
other of the eigenvalues of the matrix S^^M must ap- 
proach unity for the occurrence of a surface state or res- 
onance. This can be expressed as 

1 1 — A„ I — >■ minimum. (4) 

where A^, is the vth eigenvalue of the matrix S^^M. For 
computational purposes it is convenient to re-express this 
condition in the form 

|det [I - S"M]| minimum (5) 

where I is the identity matrix. This is the condition for 
sustained multiple scattering between the crystal atomic 
layers and the surface barrier and corresponds to the elec- 
tron being trapped in the surface region. Hence mapping 
the values of {E, k^^ ) for which Eq. (5) applies gives the 
surface band structure of the crystal surface. 

Echenique, Pendry and co-workers^'' have used Eq. (1) 
to explain the occurrence of surface image states that 
arise from the image tail of the surface barrier. Smith and 
co-workers^* and Borstel and Thorner and co-workers^^ 
have also used this condition to calculate the energy of 
Shockley and image SSs in sp-band gaps of a number of 
metals by calculating <pc from parameterised nearly-free 
electron bands. This method is not applicable gener- 
ally and because of the narrow energy range considered 
in this work, it can be expected to yield mostly trends 
in the properties of the surface barrier between different 
surfaces of the metals. 

The eigenvalues of the reduced transfer matrix R for 
the bulk crystal give the band structure of the bulk crys- 
tal in the form E{k±) for each value of fcy where k± is 
complex and is the component of the crystal momentum 
k perpendicular to a 2D atomic layer that is a plane 



of the crystal structure. This corresponds to the com- 
plex bulk band structure in the form E{k) = E(k\\, k±) 
where the usual bulk Brillouin zone (BBZ) for values of 
k is replaced by a zone in the form of a right prism. This 
prism has polygon basal planes that each are the 2D Bril- 
louin zone for fcy values in a particular crystal plane and 
sides that extend from the central plane for k± values 
of —Tr/d and n/d where d is the spacing between atomic 
layers parallel to a particular 2D crystal plane. Bulk en- 
ergy bands have finite values for Re k± with Im fcj_ = 
and have propagating wave functions. As first described 
by Heine^'^, gaps between energy bands are spanned by 
"real lines" that have "real energy lines" for Re k± and 
loops/half-loops for Im k± and have wave functions that 
are attenuating in a direction perpendicular to the par- 
ticular crystal plane since they have finite values for both 
Re k± and Im k± . Band minima that have no lower en- 
ergy maxima to which to join have real lines that extend 
to infinity. If the crystal is terminated so that the partic- 
ular crystal plane becomes the surface of the semi-infinite 
crystal the bulk wave functions with finite Im k± that at- 
tenuate towards the interior of the crystal may form the 
bulk part of a siirface state or resonance. This will oc- 
cur if the wave function amplitude and derivative from 
the bulk crystal can match the form of the wave func- 
tion through the surface and barrier potential transition 
regions and into the vacuum side of the semi-crystal. In 
the present scattering method this is achieved by detect- 
ing sustained multiple scattering in the surface region via 
the condition in Eq. (5). It is necessary to have an accu- 
rate form for the shape and height of the surface barrier 
potential including the position of the image tail and its 
saturation closer to the surface as well as realistic ion- 
core and interstitial potentials in the surface and bulk 
regions. 

The surface energy band structure £^(fc||) contains 
SS/SRs and the continuum of surface-projected bulk 
states from the £'(A;||, A;_l) band structure. In the surface 
projection of bulk bands E{k±) for values of fey, abso- 
lute, symmetry or partial bulk energy-gaps arise from the 
critical points (maxima and minima) of the bulk bands. 
Surface-projected partial bulk energy-gaps occur where 
projected bulk band gaps are also overlaid with propa- 
gating bulk bands. Surface states lie in energy ranges of 
symmetry or absolute surface-projected bulk-band gaps 
below the vacuum level, Ey. Surface resonances lie in 
ranges of surface-projected partial (quasi) bulk gaps be- 
low Ey or they have energies above E„ where they are 
degenerate with vacuum states. One very useful expres- 
sion for the wave function of an SS/SR in terms of an 
expansion in terms of bulk states that are represented 
in the linear-combination-of-atomic-orbitals approxima- 
tion, is given by BerteP^ as 
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j k± n I 



(6) 



The factor J^i^niXii''' ~ ^j) gives the atomic-orbital 
composition (s, p etc) centred at lattice sites -R^ of bulk 
band n. The +z coordinate is the inward normal to the 
surface and kj_ = n + ifi. The third factor gives the rel- 
ative contribution to the wave function from all k± val- 
ues from all bulk bands n together with the exponential 
damping factor of that contribution, Im k± = ii towards 
the interior of the crystal from the real lines and fi must 
be negative in the present scheme. A small (large) value 
of indicates a long (short) decay length for each contri- 
bution to the wave function and each contribution decays 
by a factor exp{—\fi\d) between atom surface layers with 
spacing d normal to the surface. The contribution that 
is the most slowly attenuating in space has the smallest 
value of /iniin and = l//imin is often taken as a 
measure of the decay length of the SS/SR in the bulk. 

The complex electron self-energy E expresses the 
many-body dynamic electron-electron interaction. For 
calculating bulk and surface excited states with higher 
electron energies, the E and k variation of S from ground 
state values becomes significant and must be taken into 
account. Re S alters the energy position of one-electron 
bands and Im E modifies the bands due to inelastic pro- 
cesses from single particle and surface and bulk plasmon 
excitations. These variations can be included in a phe- 
nomcnological way. Inelastic effects due to surface and 
bulk electron-phonon interaction can also be included via 
imaginary contributions to the scattering phase shifts. 

The scattering solution for the electron with (i?, fc||) 
has the wave function of the electron in the thin 2D pe- 
riodic layer of constant potential between the top atomic 
layer and the surface barrier described by forward and 
backward travelling plane waves that scatter between the 
barrier and crystal atomic layers. For the case of negligi- 
ble inelastic electron-electron interaction, the respective 
components of their wave vectors, kf parallel to the sur- 
face and along the inward surface normal of the crystal, 
for each plane wave v are given by^® 

fcj = (fc|| + f, ±fcj_,„) and (7) 

k±,v = (i^ - |fe|| (positive root, (8) 

real or imaginary) 

using Rydberg atomic units (1 Ry = 13.6 eV, 1 a.u. = 
0.5292 A). Plane waves directed inward towards the crys- 
tal have fe+ and v labels the reciprocal-net vectors of 
the surface. In the present case we have no variation 
of the potentials of surface layers of atoms with respect 
to bulk layers. The surface barrier is one-dimensional 
with only variation in the direction normal to the sur- 
face. The height of the barrier, Uq is £"„ — Ec- The 
determination of the complex reflection coefficients for 



scattering at this surface barrier involves only the energy 
associated with the motion perpendicular to the surface 
E±^v = (Re fc_L_„)^ and the wave vector normal compo- 
nents ±fc_L^„ of the plane waves. 

Historically surface states (and resonances) have been 
designated as Shockley type or Tamm type. Shockley 
type occur when hulk bands cross on forming surface- 
projected band gaps and the wave functions interact 
causing a hybridized charge density. They are sensitive to 
the surface barrier potential shape and do not require an 
unsymmetrical termination of the bulk periodic part of 
the potential at the surface. Tamm type can occur when 
there is an unsymmetrical termination of that periodic 
potential at the surface and do not require band cross- 
ing. They retain the same characteristics as the band 
from which they derived. 

In a scattering picture for this present model there are 
only two differentiated ways of forming an SS/SR. If the 
sustained multiple scattering condition of Eq. (5) involves 
some of the plane waves of the set with finite Re fcj,,^, and 
E±^„ reflecting from points along the shape of the surface 
barrier then the SS/SR main wave function amplitude 
may extend well past the centre of the top atomic layer 
of atoms into the vacuum. These SS/SRs are sensitive to 
the surface barrier shape. Also included here are Rydberg 
or image SS/SRs where reflection occurs at the image 
tail of the surface barrier. Here the SS/SR main wave 
function amplitude extends even further into the vacuum. 
We term these as surface barrier or type-one SS/SRs. 
The special case of the Shockley SS/SRs is of this type. 

If the condition of Eq. (5) involves only those plane 
waves of the set with Re k±^v = = E±^v (only finite Im 
fc_L,ii) reflecting from the bottom of the barrier then the 
SS/SR main wave function amplitude does not extend 
far beyond the top row of atoms into the vacuum. These 
particular plane waves are exponentially attenuating or 
evanescent. The bottom of the barrier corresponds in 
energy to the constant potential energy of the top atomic 
layer and spatially is where this potential just starts to 
turn upward to produce the barrier potential. Here the 
SS/SRs arise because of the unsymmetrical termination 
of the bulk periodic part of the potential at the surface 
due to the overlapping up-turn of the potential of the 
surface barrier. These SS/SRs are not sensitive to the 
surface barrier shape but only the start of the overlapping 
potential up-turn. We term these as type-two SS/SRs 
and they are more Tamm-like. 

For the scattering method here it is apparent that la- 
belling an SS/SR as either Shockley or Tamm type is 
only possible in some special cases. 

The calculation of the matrices R, Mi, M and S^^ 
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is the same as that used by McRae in the calculation of 
intensities in low-energy electron diffraction (LEED).^^ 



III. RECENT SURFACE BAND 
CALCULATIONS FOR AL(lll) 

The most complete and recent calculations of surface 
band structure for Al(lll) are an ab initio self-consistent 
finite-layer slab calculation using the density functional 
(DF) method with local-density-approximation (LDA)^ 
that was also extended to higher energies rising a non- 
local exchange-correlation (xc) potential to produce the 
image tail of the surface barrier potential®; an ab initio 
self-consistent surface-embedded Green-function (SEGF) 
method for a semi-infinite crystal^ and an ab initio self- 
consistent assembly of boundary-controlled monolayers 
(ABCM) method^" for a semi-infinite crystal. 

The DF calculations allowed for free-form surface po- 
tentials, the SEGF method used a full-potential with 
warping terms in the interstitial region and non-spherical 
terms inside the atomic muffin-tins and the ABCM 
method used the warped-muffin-tin approximation where 
there is no restriction on the shape of the self-consistent 
potential outside the muffin-tin spheres. 

From Fig. 5 of Ref. 10 we note that the valence charge 
density of the sub-surface atomic layers is the same as 
the bulk, and as noted in Ref. 10 "considerable charge 
transfer is essentially limited to the vaciium side of the 
ion cores and the charge density on the bulk side is little 
different from the bulk". The charge density shown in 
Fig. 2 of Ref. 8 is similar to that of Ref. 10. 



IV. DETERMINATION OF BULK, SURFACE 
AND BARRIER POTENTIALS. 

The lattice constant for fee Al at room temperature 
of 300 K is 4.0496 A and the (111) surface has primi- 
tive net vectors of length 2.863 A at an angle of 120°. 
The spacing between bulk layers parallel to the surface 
is d = 2.338 A. The primitive reciprocal-net unit cell is 
rhombic with sides of length 2.534 A"-'^ at 60° and the 
surface Brillouin zone (SBZ) is hexagonal. Recent LEED 
determinations of surface structure have found an expan- 
sion of the surface layer position normal to the surface 
of ~ 1.4 %^'^ and smaller variations from bulk values for 
sub-surface layers. These variations are considered to 
have a negligible effect in the present calculation and for 
clarity in the analysis of results they are not included 
here. For the same reasons electron-phonon interactions 
at 300 K are also not included. 

The ion-core scattering phase shifts were calculated 
from the bulk fee Al potential found from the self- 
consistent band structure calculation of Moruzzi et al.^'' 
This potential is of the muffin-tin form with the self- 
energy xc potential calculated from the Hedin-Lundqvist 
LDA functional. The muffin-tin radius is 1.42 A. The 



Fermi energy, Ef is 8.6 cV above the bulk constant po- 
tential energy, Ec- The symmetry designations of all 
the bands are shown in the bulk band structure calcu- 
lation of Ref. 25. More recent ab initio calculations have 
been performed for the bulk band structure of Al with 
full potentials that include non-muffin-tin corrections and 
methods to account for the self-energy S from many- 
body electron-electron interaction. The self-consistent 
DF calculation using the CRYSTAL software with use of 
a local and non-local xc functional gave results extend- 
ing to 14 eV above Ef?^ A DF calculation using the 
ABINIT software incorporated a many-body perturba- 
tion approximation to account for the self-energy. This 
GW approximation used the Green function propagator 
G and static self-consistent screening interaction W to 
give results to 12 eV above -E/.^^ For energies up to 9 
eV above Ef both of the latter calculations show small 
differences from each other and small differences with the 
band structure from Ref. 24 and that calculated in the 
present work except for the sd-likc Xf point. The above 
two DF calculations place the Xf point ~ 0.5 eV higher 
than the other two calculations. 

The surface and bulk plasmon excitation energies for 
this surface are at 10.2±0.2 eV and 15.0±0.2 eV^^ above 
Ef where inelastic electron-electron interaction becomes 
significant. Below these energies the effects of inelastic 
processes on bands from Im T,{E, k) are negligible and 
are not included in the present calculation except as a 
check in one particular case as noted later. Likewise the 
corrections necessary for excited high-energy bands from 
Re Ti{E, k) are not included here because the present 
energy range is only 9 cV above Ef. 

For simple metals^'^'^^ the reduced atom coordination 
at the surface results in band narrowing in the surface 
layer, a decrease in surface ion-core binding energies and 
decrease in the magnitude of the average potential in the 
surface layer. The spill-out of electrons at the surface re- 
sults in the opposite change in the ion-core binding ener- 
gies and magnitude of the average potential in the surface 
layer. The experimental surface core- level (binding- 
energy) shifts (SCLS) measure the initial-state contribu- 
tion from the net above effects and also a final-state con- 
tribution from the relaxation energy due to the change in 
screening between bulk and surface. For simple metals 
the final-state contribution is small^^'^^ and hence the 
experimental SCLS are approximately equal to the ini- 
tial state SCLS. Aldcn et al^'* have found that the initial 
state SCLS for simple metals is approximately equal to 
the change in magnitude of the surface-layer average po- 
tential energy. For Al(lll) the experimental SCLS of the 
2p core state for the first and second surface atom layers 
are only — 0.027±0.03 eV and respectively.^^ Hence as- 
suming the final-state SCLS to be small we approximate 
the surface ion-core potential to be the same as the bulk 
and also the interstitial potential in the surface layer to 
be the same as the bulk up to the centre of the top row of 
atoms. Past this point we add the one-dimensional sur- 
face barrier potential that arises from the spill-out of the 
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valence electrons into the vacuum region. It is to be noted 
that this effective overlapping of the spherically symmet- 
ric mufRn-tin potential at the surface and the surface 
barrier normal to the surface gives a non-spherical con- 
tribution to the surface potential. This model is also con- 
sistent with the charge density results from self-consistent 
calculations for Al(lll) using free-form surface potentials 
that were discussed above in Sec. III. 



It has an image tail of the form l/[2(2; — 2:0)] where zq is 
the position of the image plane. This joins smoothly to 
a cubic-polynomial-typc saturation at a point zi closer 
to the metal surface and also smoothly to the constant 
potential of the surface layer, i?c,i at Z2- Here we have 
set the barrier to join to the constant potential energy 
Ec^i at the centre of the atoms in the top layer at Z'2 
= 0. With the same constant potential for the surface 
layer as the bulk and energies with respect to Ec, the 
height of the surface barrier [/q = - = 12.84 cV. 
From the smooth joining of the segments and that C/q 
here is pre-determined from the bulk band structure and 
experimental work function we have only two indepen- 
dent parameters - the image plane zq and the point of 
onset of saturation Zi. The parameters zq and zi are to 
be determined from the experimental energies of SS/SRs 
that depend on the shape of the surface barrier potential 
at r. 

At r an SS has been detected by high-resolution 
ARPES at 4.56 ±0.04 eV below Ef.^^ The first Rydberg 
image resonance has been detected at 3.75 eV above Ef 
or ~ 0.5 eV below the vacuum level by k-resolved inverse 
photoemission spectroscopy (KRIPES)^^'^^ and scanning 
tunnelling spectroscopy (STS)^^ and by two-photon pho- 
toemission (2PPE)'^'' at 0.46 ±0.1 eV below the vacuum 
level. 

The shape details of the empirical surface barrier that 
are to be used in the calculations for the whole SBZ are 
found by varying the parameters zq and zi to obtain a 
match between the calculated and experimental surface 
bands at T. 



V. CALCULATION OF SURFACE 
STATES/RESONANCES FOR T. 

First we calculate the complex bulk band structure in 
the form E{k) = ^^(feii, fcj,) where fey = corresponds to 
r of the SBZ from the bulk potential from Ref. 24. This 
complex bulk band structure also shows attenuating bulk 



The experimental work function $ for Al(lll) is 4.24± 
0.02 cV^^ and this gives the vacuum level E^ with respect 
to Ec&s E^-Ec^ Ef + <^ = 12.84 eV. We use an empir- 
ical surface barrier shape Ui){z) that has been suggested 
by Rundgren and Malmstrom^*^ where z is the perpen- 
dicular distance from the centre of the first row of atoms 
which is located at 2 = and the barrier is located at 
negative z values. In Rydberg atomic units it is 



(9) 

I 

states. The energy interval in all the calculations in this 
work is 0.01 eV. Seven plane waves were found to be ad- 
equate for the plane- wave basis set for interlayer scatter- 
ing at r. Only the plane wave with reciprocal-net-vector 
coefficients 00 is propagating in this energy range and 
the 3-fold degenerate {01} and {10} sets are evanescent 
plane waves. Spherical partial waves up to I max = 5 
were sufficient to describe the intralayer scattering. The 
resulting band structure is shown in Fig. 1(a) and (b) and 
passes through the F and L points of the BBZ as shown. 
All energies are subsequently referred to Ej unless stated 
otherwise. Hence we see that the gap between bulk bands 
(blue/grey) at L2' and Li near —4.6 eV projects as an 
absolute surface gap at F. We also note the real lines with 
real energy lines (orange/light grey) in Fig. 1(a) in the 
sp L gap and small loop (orange/light grey) in Fig. 1(b) 
that is also shown in the insert on an expanded scale. 
Bands (blue/grey) in Fig. 1(a) give lines at Im ~ Q 
only (blue/grey) in Fig. 1(b). We also note the real en- 
ergy lines and half-loops that start from bulk minima at 
higher energies and extend to —00. These are from min- 
ima in A3 and Ai bands that are shown in Ref. 25. These 
real energy lines and half-loops indicate a very wide pro- 
jected partial even-symmetry gap from the Ai minimum 
(15.5 eV) to —00 and a wide odd-symmetry gap from the 
A3 minimum (10.5 eV) to —00. Other real energy lines 
and half-loops with larger Im k± values for the present 
energy range from higher energy minima are not calcu- 
lated because they have negligible effect at these lower 
energies and also lead to instabilities in the calculation. 

Using Eq. (5) we calculate the energies of the SS/SRs 
at r where fcy =0 and find that the first Rydberg im- 
age SR energy is reproduced at —0.50 eV with respect to 
Ey with Zq at —0.53 A and z\ at —1.06 A. Also without 
any change in these Zq and Zi values the lower SS is at 
—4.58 eV and in the correct experimental energy range. 
At the condition of Eq. (5), only the 00 plane wave is 
propagating with finite Re k±_ and is incident at points 
along the shape of the surface barrier and reflects from it. 
Hence these SS/SRs are type one. The SS conforms to 



Ub{z) = 




, + l/[2(z-zo)] 

+ ^1(2; - Zi) + V2{Z - Zif + Vs{z - Zif 



for Z < Zi < Zq, 

for zi < z < Z2, 
for z > Z2- 
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FIG. 1: (color online) The real and imaginary parts of the 
complex bulk and surface energy band structure -E(fc||, k±) at 
fc|| = = r for Al(lll) are shown in panels (a) and (b) respec- 
tively. These k values pass along a line through the F and L 
points of the BBZ as indicated in panel (a). Propagating bulk 
bands are blue / grey points and real lines of attenuating bulk 
bands are orange/light grey points. The real lines that extend 
to —CO from 'unmatched' bulk-band minima are terminated 
in the diagram at low energies where they become unstable 
and have negligible effect in the calculation. Other real lines 
from such minima at higher energies that do not produce a 
significant effect in the calculation and have large Im k± val- 
ues are not shown. Bulk band labels at L are indicated at the 
margins. The insert in panel (b) shows an expanded scale for 
the energy range —5 to —4.2 eV. Note that loops/half-loops 
in panel (b) show most readily the surface-projected absolute 
and partial bulk band gaps. Vertical full black lines indicate 
energies of surface states or resonances. The intersection of 
vertical dashed black lines in panel (b) with bulk Im k± values 
indicates attenuation factors \fi\ of each bulk-band contribu- 
tion to the surface state or resonance. Near 4 eV are the 
first, second and last (n = oo) members of the Rydberg se- 
ries of image surface resonances. Panel (c) shows the result 
of Eq. (5) for the above fc|| value where minima indicate en- 
ergies of surface states and resonances. The minium at —8.6 
eV corresponds to a surface resonance for a free electron in 
the crystal. 
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FIG. 2: (color online) Surface barrier potential energy rela- 
tive to Ec, the bulk constant potential energy for Al(lll). 
The dashed line is the continuation of the unsaturated image 
potential. The image plane position zo — —0.53 A and satu- 
ration begins at zi = —1.06 A. The muffin-tin radius is 1.42 
A and the geometric surface Zj = —1.17 A. The vacuum and 
Fermi level are Ev and Ef respectively relative to Ec. 

the classical Shockley type. It is placed slightly closer to 
the Li gap edge in the L2' , Li absolute gap with width 
calculated here to be 0.29 eV. The position of the SS 
with respect to gap edges is not certain experimentally 
and there are no experimental determinations of the en- 
ergy positions of the L2', Li points. As the accuracy 
of the calculated positions of the L2', Li points from the 
present band structure cannot be assessed we have not 
attempted to also position the SS with respect to either 
gap edge. The result from Eq. (5) is shown in Fig. 1(c). 
This result was found previously where it was used to de- 
termine SRs above the vacuum level to 27 eV at F where 
inelastic scattering from Im E is significant and must be 
included. Here we are interested in extending the calcu- 
lation to all fc|| points of the SBZ but for energies where 
Im E is negligible. Scattering phase shifts at the image 
part of the barrier in multiples of 27r from that found for 
the above zq will position the image SRs at almost the 
same energy. The values zq ~ —2.70 A and zi = —4.2 A 
also produce the correct first Rydberg image SR energy 
and also the correct SS energy. However these values also 
produce an SR at 0.60 eV at F for this wider barrier. At 
present there is no experimental evidence for such an SR. 
Hence the first values for the shape of the surface barrier 
were used in the present calculation. This barrier shape 
is plotted in Fig. 2. 

The above calculated energy positions of the SS and 
SR, the first, second and last (n = 00) Rydberg image 
SRs and projected band gap for F are plotted on the sur- 
face band structure E{k\\ ) for M - K - T - M shown in 
Fig. 3. The SS/SRs are also plotted in Fig. 1(a) (vertical 
black lines) where they have no dispersion with values 




2.1 2 1.9 1.8 1.7 1.6 1.5 1.4 1.2 1 0.8 0.6 0.4 0.2 0.2 0.4 0.6 0.8 1 1.2 

^„(A-') fc„(A-') ^„(A"') 

FIG. 3: (color online) Surface energy band structure i5(fe||) for Al(lll). Surface states and resonances are indicated by full 
black lines. Thick full lines indicate surface states or (sections of) strong surface resonances where substrate reflection is strong. 
Thinner full lines indicate (sections of) weaker surface resonances. The continuum of surface-projected bulk bands are shaded 
blue/grey and their absolute gaps unshaded. Dashed black lines indicate the free-electron surface band structure labeled with 
the coefficients of the surface reciprocal-net vectors. All SS/SRs below the lowest energy free-electron band are necessarily type 
two. The horizontal axes are scaled differently in order to display the relevant features more clearly. 



of Re k± , to show the full surface and bulk band struc- 
ture at r. These full band structures for a given value of 
fc|| show surface, propagating bulk and also attenuating 
bulk bands near the surface and all are important in the 
analysis of energy and momentum distribution curves in 
photoemission and in other surface spectroscopies. The 
attenuation constants = |Im k±\ of each bulk-band 
contribution to the SS/SRs are indicated from the SS/SR 
energy intersection (vertical dashed black line) with bulk 
Im k± values. The small values of Im k± in the loop 
in Fig. 1(b) for the L gap and at the energy of the SS 
account for the fact that this SS has a wave function 
component with large decay length and long penetra- 
tion into the bulk of the crystal as well as components 
that do not penetrate much past the surface layer. From 
Fig. 1(b) As = l//imin ~ 71 A in the present calcula- 
tion. This SS is related to the SS found experimentally 
on Al(lOO) at -4.55 eV at X in the absolute L gap.^S'^i 
For the bulk and surface potentials used here the SS lies 



slightly closer to the upper Li (sii-like^^) band edge at 
-4.49 eV than L2' (p-like^^) at -4.78 eV and there is an 
SR just above the Li gap edge at —4.45 eV. If the surface 
barrier parameters are changed slightly so that the SS is 
formed nearer the lower energy gap-edge then this extra 
SR changes position to just below the L2' edge. This is a 
very surface localised SR with calculated As ~ 0.5 A that 
may affect photoemission energy and momentum distri- 
bution curves. This SR is in the very wide partial pro- 
jected band-gap formed between the band minima from 
the higher energy Ai bulk band and —00 that is men- 
tioned above. The Rydberg image SRs are also in this 
wide partial gap with calculated As ~ 0.7 A. This type of 
partial projected gap will be discussed again later. The 
separation between the surface and sub-surface atomic 
layers here is 2.34 A. 
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FIG. 4: (color online) The real and imaginary parts of the 
complex bulk and surface energy band structure -E(fc||, k±) at 
fc|| = M for AI(lll) are shown in panels (a) and (b) respec- 
tively. These k values pass along a line through the L and 
X points of the BBZ as indicated in panel (a). Bulk band 
labels at L and X are indicated at the margins. A very small 
gap at 5.3 eV between the band maximum at X^i and the 
band minimum at X3 is not discernable on the scale of the 
diagram. Other features of panels (a) and (b) are the same 
as Fig. 1. Panel (c) shows the result of Eq. (5) for the above 
fc|| value where minima indicate energies of surface states and 
resonances. The insert shows an expansion for the energy 
range 1 to 1.3 eV. The minium at —2.49 eV corresponds to a 
surface resonance for a free electron in the crystal. 



VI. CALCULATION OF SURFACE 
STATES/RESONANCES AT M. 

The calculated complex bulk band structure for 
E{ki\,kj_) with = 1.27 A^i EE M/m'(M{L)/M{X)) 
point is shown in Fig. 4(a) and (b) that also shows 
attenuating bulk states. Note that loops/half- loops in 
panel (b) show most readily the surface-projected abso- 
lute and partial bulk band gaps. We find convergence 
with up to eight plane waves in our region of constant 
potential between scattering layers. Here the bulk bands 
pass through the L and X points of the BBZ as shown. 
The L2', Li bulk gap near —4.6 eV becomes a surface- 



projected partial bulk band gap at M because it is over- 
laid with other free-electron-like bulk bands in the sur- 
face projection. The bulk gap between the X41 point at 
—2.98 eV and Xf point at —1.98 eV is also a projected 
partial gap overlaid with free-electron-like bulk bands. 
Real lines spanning these bulk gaps are real energy lines 
in Fig. 4(a) and loops in Fig. 4(b). Energy bands in 
Fig. 4(a) give lines at Im k± — only in Fig. 4(b) that 
also pass through loops in the case of projected partial 
gaps. The 2-fold degenerate X^r point is at 5.30 eV and 
the X3 point is slightly higher with a very small 0.02 eV 
gap between these points giving a projected partial gap 
between 5.30 and 5.32 eV. The real energy line and loop 
for this partial gap are not discernable on the scale of 
these diagrams. The X" point is at 8.40 eV. From X^i 
in Fig. 4(a) the band that has a maximum there and dis- 
perses downward has even bulk symmetry and the other 
flatter band that is a minimum there and disperses up- 
wards has odd bulk symmetry. The band minima at X^ 
(8.40 eV) and at X5, (5.30 eV) from the flatter band 
have no maxima to which to join and they produce real 
energy lines for Re k± and half-loops for Im k± extending 
to —00. Since all bands below X5' have even bulk sym- 
metry, a wide projected odd-symmetry bulk gap exists 
from X51 to —00 and a wide even-symmetry partial bulk 
gap from X" to —00. An absolute projected gap exists 
near 1 eV and below the band minimum at —5.37 eV. 

The result of the calculation for the SS/SRs from 
Eq. (5) is shown in Fig. 4(c) and energies plotted in 
Fig. 4(a) & (b) and in Fig. 3 for M together with the 
absolute projected gaps. The minimum at —2.49 eV in 
Fig. 4(c), as well as the minimum at —8.6 eV at F in 
Fig. 1(c), are not SRs of this system but indicate the crys- 
tal empty-net or surface free-electron energy. This occurs 
when a plane wave in the region of constant potential be- 
tween layers transitions from evanescent to propagating. 
These energies are also plotted in Fig. 3. Above —2.49 
eV only the two degenerate plane waves with reciprocal- 
net-vector coefficients 00 and 10 become propagating in 
this energy range. Below —2.49 eV these plane waves and 
all others are evanescent and have Re fc_L = and finite 
Im k±. The results in Fig. 4(c) determine that there are 
three occupied SRs below —2.49 eV at M. Two are at 
—4.80 and —4.48 eV which are near the energies of the 
partial gap near —4.6 eV with edges at I/2' (—4.78 eV) 
and Li (—4.49 eV). Both of these SRs are just outside 
the edges of this partial gap and have small As values 
« 0.6 A. However a small change in the partial gap po- 
sition and width from the bulk band structure could put 
one or both of these SRs inside the partial gap and con- 
sequently with a large As value. Another SR is at —2.96 
eV in the partial gap with edges X4/ at —2.98 eV and Xf 
at —1.98 and 0.02 eV above the lower energy p-like X41 
edge. The calculated As is « 53 A. This SR is related to 
the SS observed on Al(lOO) at -2.75 eV at F in the ab- 
solute X gap and ~ 0.1 eV above the X41 energy.^*'^^''^^ 
The three SRs here arise from scattering at the bottom of 
the surface barrier and are not sensitive to its shape but 
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are slightly sensitive to its starting position Z2- Accord- 
ing to our classification in Sec. II they are type-two SRs 
with wave function maximum amplitude not extending 
far beyond the top row of atoms into the vacuum. 

There is no identification of occupied SRs at M /m' in 
the ARPES experimental data of Hofmann and Kambe^^ 
and no other experimental data at M /M . 

Above —2.49 eV the 00 and 10 plane waves have finite 
Re fc_L and are incident at points along the shape of the 
surface barrier potential and reflect from it. The absolute 
gap near 1 eV has edges at 0.80 and 1.26 eV and has two 
SSs at 1.14 eV (A^ « 41 A) and 1.17 eV (A^ « 48 A) 
and there is an extra SR just above this gap with small 
Ag. This SR is in the wide partial gap between the band 
minimum at X" and — oo. There is also an SR at 5.30 
eV in the very small X5/ , Xj, partial gap with large A^ . 

A KRIPES experimental study of unoccupied surface 
bands in this region found an SR at 4.0 eV for M 
whereas we calculate an SR at 5.30 eV. Hence there is a 
1.3 eV discrepancy. M and Af are different only with re- 
spect to possible allowed initial and final state transitions 
in (inverse) photoemission. The bulk band structure in 
Figs. 4(a) and (b) for M differs from that of M only that 
the sign of k±^ values are interchanged. The partial gap 
at M and M near 5.3 eV that we calculate is at the same 
energy as the partial gap in Fig. 6(b) at M in Ref. 37. 
We also find the dispersion of this gap for F - m' to be 
the same but Fig. 6(b) of Ref. 37 shows a projected gap 
for M at ~ 1.8 eV whereas ours is centred at 1.0 eV. The 
Al band calculations mentioned previously that use cor- 
rections for electron self-energy S for excited states^^'^^ 
both place the X^i point near 5.5 eV although the gap 
to ^3 is wider than that found here. Inelastic effects in 
our calculation are considered to be insignificant because 
the energy is well below the onset of surface plasmons 
at 10.2 eV and bulk plasmons at 15.0 eV. However as 
a check we can include these effects from Im Yi{E,k) in 
our calculation according to Eq. (1) and Fig. 4 of Ref. 35. 
We find that this inclusion has negligible effect on the en- 
ergy position of this SR. The other barrier model with 
= —2.70 A and zi = —4.2 A also produces the SR 
at this same energy of 5.3 eV. The surface band calcula- 
tion of Ref. 10 extending to 5 eV above Ef also did not 

produce an SR at 4.0 eV at M/m'. 



VII. CALCULATION OF^URFACE BANDS 
FOR r - M. 

The above methods were used to calculate surface 
bands E{k\\) and the continuum of surface-projected bulk 
bands and their gaps for fcy values corresponding to F - 
M along the S direction and the results are plotted in 
Fig. 3. Approximately 40 values of fcy were used to span 
this range of fcy. Thick full lines indicate surface states 
or strong resonances. Thinner full lines indicate weaker 



surface resonances. Resonances are classified as strong if 
in the present scattering picture they are formed where 
substrate reflection is strong. The dashed lines show the 
crystal empty-net or free-electron surface band structure 
that is also the (S, fcy) values at which plane waves in 
the region of constant potential between layers transi- 
tion from evanescent to propagating in the crystal (i.e. 
acquire real fcj_). To understand the dispersion and na- 
ture of the SS/SRs the full complex band structure in 
the near-surface region for a number of key values of fcy 
for F - M is examined. 

Fig. 5 shows our full bulk and surface band structure 
for = 1.21 A~^. Particularly notable is the gap near 1 
eV at M in Fig. 4 becoming two partial gaps near 1.7 eV 
and 0.5 eV with SRs following inside these partial gaps. 
Note that there are SRs with small As also following just 
outside these partial gaps. These latter SRs are in the 
wide partial gaps formed from the higher energy band 
minima at 8.4 and 5.32 eV and —00. Below Ef a new 
SR appears at —5.28 eV that is not a continuation of the 
SR at -4.68 eV. The SR at -4.68 eV and also the SR 
at —4.39 eV, have small As ~ 0.5 A and are in the same 
wide gap mentioned above while the —5.28 eV SR is in 
the less-wide partial gap between the band minimum at 
-4.8 eV and -00 with u,ni„ ~ 0.38 A"^ and As ~ 2.6 A. 

Fig. 6 shows our full band calculation for fc|| = 

1.04 A" 1 along f -M' . In their KRIPES study, Yang 
et aP^ interpret their experimental feature at ~ 4.2 eV 
to be an SR while we calculate unoccupied SRs at 3.31, 
3.43, 5.55 and 5.88 eV. The SRs at 3.43 and 5.88 eV 
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FIG. 5: (color online) The real and imaginary parts of the 
complex bulk and surface energy band structure -E(fe|| ,k±) at 
fc|| = 1.21 A"^ along r - M for Al(lll) are shown. A very 
small gap occurs at 5.3 eV that is not discernable on the scale 
of the diagram. Other features are the same as in Fig. 1(a) 
and (b). 
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FIG. 6: (color online) The real and imaginary parts of the 
complex bulk and surface energy band structure £(fc||, k±) at 
fcy = 1.04 A"^ along T -M' for Al(lll) are shown. A very 
small gap occurs at 5.5 eV that is not discernable on the scale 
of the diagram. Other features are the same as in Fig. 1(a) 
and (b). 



have small values of and are in the wide partial gap 
extending from the band minimum at 8.6 eV to — oo. 
It is worth considering that interpretations from differ- 
ences in the calculated bulk band structures could be 
involved. Reference 37 shows their bulk band calculation 
for this fc|[ in their Fig. 6(a). They use an /-dependent 
self-consistent pseudo-potential method that is different 
from our method. However the two bulk band structures 
are nearly identical so that this is not a factor in the dis- 
crepancy. Note that there are image SRs near 8 eV also 
in the above wide partial gap with As values of « 5.3 A 
that are much larger than those at P. Below Ef the two 
lower energy SRs at —2.37 and —1.29 eV are also in this 
wide partial gap with small As while the SR at —0.60 eV 
is in a small width partial gap with As ~ 35 A. 

Fig. 7 shows our full band structure for fc|| = 
3/4(rM) = 0.95 A-i. The bulk bands pass through 
the K and U points of the BBZ. The gap between the 
minima and maxima of the bands near 5.78 eV is still not 
visible on the scale of this diagram but the flatter odd 
symmetry band can be seen to have its minimum moved 
to 5.64 eV. There is an SR at 5.78 eV in the very small 
gap and an SR at 6.10 eV with As ~ 1.7 A in the larger 
gap extending from the band minimum at 8.9 eV to —oo. 
Two partial gaps of small width below Ef in Fig. 6 now 
become one partial gap containing one SR at —1.90 eV 
just outside this partial gap with small As and one at 
-0.86 eV (As ~ 21 A) inside it. 

In Fig. 8 for fc|| = 0.92 A~^ a new very small par- 
tial gap appears at —1.8 eV that moves down in energy 
with decreasing ku and persists to become the absolute 
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FIG. 7: (color online) The real and imaginary parts of the 
complex bulk and surface energy band structure E{k^^,k±) 
at fc|| = 0.95 A"^ = 0.75(rM) along F - If for Al(lll) are 
shown. These k values pass along a line through the K and U 
points of the BBZ as indicated in panel (a). Bulk band labels 
at K and U are indicated at the margins. A very small gap 
occurs at 5.78 eV that is not discernable on the scale of the 
diagram. Other features are the same as in Fig. 1(a) and (b). 



L gap at r in Fig. 1 and the SR from Fig. 7 at -1.90 
eV continues as an SS in this case. The other partial 
gap with small width below Ef has moved up in energy 
from its position in Fig. 7 along with the SR in it. At 
< 0.75 A"-'^ for r - M this partial gap disappears and 
so does the SR that has followed it as shown in Fig. 3. 

There have been a number of studies of SS/SRs below 
Ef for r - M from ARPES. Kevan et api found the SS at 
P dispersed upwards and terminated at —2.0 eV for fc|| w 
0.9 A^^ and this agrees with our result shown in Fig. 3. 
Apparently they did not detect the SS continuation as an 
SR. Grepstad and Slagvold^* measured an SR dispersing 
from _-0.4 eV at 1.08 A'^ to -0.7 eV at 0.95 A'^ near 
3/4(PM) that agrees with our result. This SR dispersion 
was also measured by Hofmann and Kambc.^^ In addition 
the latter find what they refer to as "partial gap-edge 
features" from -3 to -2 eV at k\\ = 1.15 A^^ to 0.90 A^^ 

along P - m' and also within this fey range near —1.2 eV. 
They do not detect these features along P - M and hence 
do not consider them to be SRs. They consider them to 
be "a kind of surface enhancement occurring at partial 
gap edges". However they coincide with our SRs in this 
range that have very small As ~ 0.7 A and are in the 
wide partial gap from the higher energy band minima 
and —00 but just outside a smaller partial gap as shown 
in Figs. 6 and 7. These two SRs are very surface localized 
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FIG. 8: (color online) The real and imaginary parts of the 
complex bulk and surface energy band structure -E(fc||, k±) at 
fc|| = 0.95 A"^ along T - M for Al(lll) are shown. Other 
features are the same as in Fig. 1(a) and (b). 



and appear to be difficult to detect in some theoretical 
methods and will be referred to later. 
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VIII. CALCULATION OF SURFACE 
STATES/RESONANCES AT K 

The calculated complex band structure for ^^(fciijfcj^) 
and fc|| = 1.46 = K point is shown in Fig. 9(a) 

and (b) that also shows attenuating bulk states. Here 
there are five absolute projected gaps; two above and 
three below Ef including the absolute gap below the band 
minimum at —3.5 eV. Real energy lines and half-loops 
extend from a band minima ~ 15 eV and produce a wide 
partial projected gap from this energy to — oo. The result 
of the calculation for SS/SRs from Eq. (5) is shown in 
Fig. 9(c) and energies plotted in Fig. 9(a) & (b) and in 
Fig. 3 for K together with the absolute projected gaps. 
The minimum at —0.45 eV in Fig. 9(c) is not an SS of 
this system but indicates the crystal empty-net or surface 
free-electron energy as explained earlier. Above —0.45 eV 
only the three degenerate plane waves with reciprocal- 
net-vector coefhcients 00, 10 and 01 become propagating 
in this energy range. Below —0.45 eV these plane waves 
and all others are evanescent and have Re k± = and 
finite Im fc_L. The results in Fig. 9(c) determine that there 
are three SSs and two SRs below Ej that are all type two 
in our classification. The SSs in the lowest narrow gap 
are at -2.87 and -2.60 eV with of f« 23 A and 53 A 
respectively. The SS at -0.70 eV has « 11 A. The 
SRs are not far outside the absolute gap edges with very 
small As w 0.4 A in the very wide partial gap from sa 15 



FIG. 9: (color online) The real and imaginary parts of the 
complex bulk and surface energy band structure -E(fe|| ,k±) at 
fc|| = K for Al(lll) are shown in panels (a) and (b) respec- 
tively. Other features of panels (a) and (b) are the same as 
Fig. 1. Panel (c) shows the result of Eq. (5) for the above 
fc|| value where minima indicate energies of surface states and 
resonances. The insert shows an expansion for the energy 
range 3 to 3.3 eV. The minium at —0.45 eV corresponds to a 
surface state for a free electron in the crystal. 

eV to — oo. Above —0.45 eV the three degenerate plane 
waves labelled 00, 10 and 01 have finite Re k± and are 
incident at points along the shape of the surface barrier. 
Above Ef there are three SSs at 3.09, 3.13 and 3.20 eV 
with As ~ 31, 33, 38 A respectively in one gap and two 
at 7.02 and 8.31 eV with As « 66 and 16 A respectively 
in the higher energy gap. Note that there are also two 
SRs that are close to the absolute gap edges with very 
small As ~ 0.45 A again in the partial gap from 15 eV 
to —oo. 

Only the SS at —0.7 eV has been fomid from the 
ARPES experiments of Hofmann and Kambe^^ and Ke- 
van et al.^^ 



IX. CALCULATIONS FOR T - K - M 

Approximately forty fcy values corresponding toT - K 
along the T direction were used to produce the surface 
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band structure plotted in Fig. 3 from the calculated com- 
plex bulk band structures and using Eq. (5). As this is 
not a bulk symmetry direction there are many SR bands 
degenerate with the bulk bands and the surface band 
structure is complicated. Kevan et al^^ found the —0.7 
eV SS at K disperses upward towards E j for K - T over 
a range of fc|| « 0.3 A^^ and our calculated dispersion 
agrees with this. If one continues fcy values along the F 
- K direction then one arrives at M in the second SBZ. 
This was done to produce the surface band structures 
corresponding to K - M along the T direction that is 
plotted in Fig. 3 using approximately 40 values of fcy and 
the same methods as before. All SS/SRs that have en- 
ergy below the lowest energy surface free-electron band 
shown dashed in Fig. 3 are type two. 

X. COMPARISON WITH RECENT AB INITIO 
SURFACE BARRIER POTENTIAL 
CALCULATIONS 

The position of the geometric surface is defined as the 
plane where the bulk crystal is divided in half to create 
the semi-crystal. This is also the surface edge, Zj for 
the jellium model of a semi-crystal. For Al(lll) this 
plane extends a distance 2.21 a.u. (1.17 A) from the 
centre of the top row of atoms at z = 0. The surface 
barrier height Uq depends on the work function and the 
Fermi level, Ef with respect to the constant potential 
energy in the surface atomic layer, Ec,i or layer mufBn- 
tin constant. Here Ec,i is the same as the bulk constant 
potential, Ec. Both Ec and i?c,i have a contribution 
from the electron self-energy S via the xc-potential in 
the crystal or atomic layer. Although we have used the 
experimental work function $ of 4.24 eV, there is still 
some uncertainty in the bulk band structure calculation 
of Ec and Ef. We have found that small differences in the 
values of Ec, Ef and Ec^i and hence Uq can significantly 
change the determination of the image plane position, zq. 
Coupled with the measured uncertainty of ±0.1 eV in the 
first Rydberg image SR position at 0.46 eV from Ey and 
no experimental measurement of the second image SR 
or position of the L2' and Li bands there is scope for 
some uncertainty in our value of Zq. As only Rydberg 
SS/SRs are very sensitive to Uq and zq this uncertainty 
in Zq has only a small effect on the energy position of 
other type-one SS/SRs. 

White et aP^ have calculated the electronic self-energy 
S experienced by electrons of the crystal in the near- 
surface region of Al(lll) in a repeated-slab geometry of 
five Al layers and eight vacuum layers. They use the 
ab initio GW method with the random-phase approx- 
imation for a dynamical screening interaction, W and 
a non-interacting approximation for the Greens func- 
tion, G from an initial density functional calculation 
with local-density approximation. Their non-local xc- 
potential gives the image potential and image-plane posi- 
tion, Zq. These authors find |2;o| = 0.4 ±0.2 a.u. from the 



geometric surface or 2.61 ±0.2 a.u. from z = with satu- 
ration starting at \zi\ « 6.0 + 2.21 a.u. from z = 0. From 
their figures one determines that they found Uq ~ 9.5 
eV and Ef ss 4.8 eV from Ey and hence $ ss 4.8 eV. 
In the following all SS/SRs are given at F with respect 
to Ey. They calculate an SS at 1.66 eV from Ey. An- 
other repeated-slab calculation by Heinrichsmeier et al^ 
also using the GW approximation and a one-dimensional 
potential with several hundred layers found \zq\ = 0.38 
a.u. from zj or 2.59 a.u. from z = and $ = 4.82 eV and 
found an image SR at w 0.4 eV. These two calculations 
of $ have a difference of ss 0.6 eV compared with the 
experimental value 4.24 eV. 

Hence the present result of |^o| = 1-1 a.u. (0.58 A) is 
less than half the GW ab initio result of 2.6 a.u. (1.4 A) 
from z — and inside the geometric surface in our case. 
This difference may be accounted for in part because of 
the differences in the barrier heights of these works com- 
pared with our value of [/q = 12.84 eV. If we set Uq = 9.5 
eV, \zo\ = 2.61 a.u. with \zi \ = 8.21 a.u. both from z = 
in our calculation we calculate SRs at 1.7 eV and 0.2 eV. 
This is in agreement with the result in Ref. 39. With the 
same zq and z\ values but with Uq = 12.84 eV we calcu- 
late SRs at 2.6 eV and 0.3 eV. There is no experimental 
detection of SRs between 1.7 and 2.6 eV. There is discrep- 
ancy between zq (and Uq) from ab initio and the present 
empirical model calculations. The value of Zq found from 
either calculation cannot be considered definitive at this 
stage. 



XI. CONCLUSION 

In comparing the results of the present and three 

rcccnt®"^'^ surface band calculations one is faced with 
explaining why these differ. All four agree on the oc- 
currence of surface states in surface projected bulk gaps. 
The ABCM^° method finds only surface resonances that 
are the continuations of these surface states into par- 
tial projected bulk-gaps that have a small energy width 
and have real lines between bulk band minima and max- 
ima. These resonances all have long decay lengths. The 
present scattering method produces resonance bands in 
close agreement with all of these bands but also produces 
additional bands. The slab DF^ and SEGF^ methods 
both do not produce these same resonance continuations 
of the states in all cases. The slab DF method docs show 
broad resonances in some of these cases. Both the present 
scattering method and the slab DF calculation find a res- 
onance band about 0.4 eV above the lower valence band 
edge along the T direction {M K) and extending also 
partially along S near M that is not found in either of the 
other calculations. In the present calculation this band 
and the one above it occur in the very wide partial bulk 
gap that contains attenuating bulk bands from real lines 
extending from a higher-energy bulk band minimum to 
—00. Both SRs are very close to the bulk band edges of a 
much smaller-width partial gap within this wider partial 
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gap as shown at M in Fig. 4. In our case of a continuum 
of bulk bands, the energies of these SRs do not also fall 
in the region of the smaller-width partial gap. However 
the finite layer slab DF method does not have a contin- 
uum of bulk bands and a small partial gap may have 
appeared in this energy region that may have elevated 
one of these bands to one with a long decay length and 
hence detectable in that method. 

It is suggested here that the above bands along M 
- K are two members of a class of bands not detected 
in the ABCM method. Attenuating bulk states exist in 
very wide partial gaps where real lines extend from bulk 
band minima to — oo when there is no band maximum to 
which to join. We have found in the present work that 
SRs can form also from the real lines of these bulk band 
minima if their wave function matches that through the 
surface potential transition region into the vacuum just 
as for other slower attenuating bulk states from small- 
width gaps. The lower the energy of the SR from the 
bulk band minimum the shorter the decay length pro- 
vided the SR energy does not also lie in a smaller-width 
absolute or partial gap. An SR in this type of wide gap 
does not always have a very short decay length and an 
example is shown for /ey = 1.04 for F - M in Fig. 6 
for the image SRs near 8 eV. Where these SRs have very 
short decay lengths their effect on photoemission inten- 
sity profiles and whether they can be detected experi- 
mentally remains to be determined. In the cases where 
these SRs are found theoretically to lie also just outside 
the edges of small energy-width partial bulk gaps it is 
possible that a small change to the bulk or surface po- 
tentials could move them inside this gap and change the 
decay length from short to long. 

It is also suggested here that unexplained features in 
experimental photoemission previously described as sur- 
face enhancement occurring at partial bulk gap-edges or 
enhancement of bulk density-of-states near the surface 
are actually SRs of the above type. This was discussed 
in Sec. VII. It is important for theoretical surface band 



structure methods to be capable of identifying all SRs 
including this type which may be very localized to the 
surface atomic layer and if type two (scattering from the 
bottom of the surface barrier potential) may not extend 
far into the vacuum either. Rydberg image resonances 
may also have very short decay lengths in the crystal as 
in the case here at F but they are type one (scattering 
near the top of the surface barrier) and they extend far 
into the vacuum. 

The present scattering method may have some advan- 
tages over other surface band structure methods. SufS- 
cient accuracy is obtained before the onset of numerical 
instabilities because of the rapid convergence obtained 
in the layer scattering calculation from direct and recip- 
rocal net summations. The locations of SS/SR bands 
depend strongly on the accuracy of the bulk potential. 
In the present method a bulk potential is chosen a pri- 
ori from a bulk band structure that has verified accuracy 
from comparison with experiment. The empirically con- 
structed surface potential has proven to be satisfactory 
since it reproduces all of the surface band structure ob- 
tained in the ah initio self-consistent warped muffin-tin 
surface potential used in the ABCM method. The scat- 
tering method identifies many SRs that have very short 
decay lengths and are mainly localized within the surface 
layer. These SRs have not been detected in the ABCM 
or other calculations. In the scattering method these SRs 
are just as easily detected as the SRs of long decay length. 
In addition it produces the image SS/SRs and it can be 
extended to unbound SRs in higher energy regions where 
inelastic effects become significant. The mcithod also pro- 
vides a single picture for understanding the formation of 
what have historically been termed Shockley, Tamm and 
Rydberg SS/SRs. It therefore provides a complete sur- 
face band structure that can aid in the analysis of exper- 
imental features from band transitions in photoemission 
and other surface spectroscopies. The use of the method 
for other metallic surfaces and those with adlayers of for- 
eign atoms is to be determined. 
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